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Superfluid-to-superfluid phase transitions in a Fermi-liquid leading to the 
emergence of two-gap superfluid states have been studied. There are con- 
sidered the systems of fermions of one and two (nuclear matter) sorts. The 
self- consistent equations describing new two-gap superfluid states have been 
obtained. These equations differ essentially from the equations of the BCS 
theory. The solutions of the self- consistent equations corresponding to one- 
gap and two-gap superfluid states have been found. 

PACS numbers: 71. 10. Ay, 71.10.Li. 

1. INTRODUCTION 

As a rule, superfluid states in a Fermi-liquid (FL) are considered as 
arising as a result of a phase transition from the normal state. However, 
phase transitions that might also take place in a superfluid FL can lead to 
the emergence of new superfluid states. Thus, we are speaking of phase 
transitions from one superfluid state to another. A new superfluid state can 
be characterized by not one but a few order parameters. In this report we 
consider the systems of fermions of one and two sorts. In the first case the 
possibility of a phase transition to the state for which the spin of a Cooper 
pair has no definite value but is equal with certain probability to zero or 
unity is studied (singlet-triplet pairing of fermions). In the second case a 
phase transition in superfluid nuclear matter to the state corresponding to 
the superposition of states with singlet-triplet and triplet-singlet pairing of 
nucleons (in spin and isotopic spaces) is considered. 
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2. SYSTEMS OF FERMIONS OF ONE SORT 

Our analysis is based on the theory of a superfluid FL 0. A Fermi 
superfluid is described by the normal f KlK2 = Tr ga+ 2 a Kl and anomalous 
9k!k 2 = Tr^o K2 a Kl distribution functions (g is the density matrix of the 
system, a+ and a K are the creation and annihilation operators of a fermion 
in the state with momentum p, spin projection a; for a nucleon in nuclear 
matter the state is characterized also by the isotopic spin projection r). 
The energy of the system, E = E(f,g), determines the fermion one-particle 
energy and the matrix order parameter 

dE _ dE 



The distribution functions / and g are related with the quantities e and A 
by the formulas i: 

/ = Ifl_i tanh MV 5 = _J_ tanh M. A; (2) 
J 2\ E 2 / ' y 2E 2 ' v; 

E = yJ? + AA+, t = e + y Q 

Here Yq = 1/T, —Y^/Yq = [i, T is temperature, fx is chemical potential. 

For a superfluid FL with the singlet-triplet (ST) pairing of fermions the 
matrix order parameter has the form 

A(p) = A (p)<j 2 + A(p)«57T 2 , (3) 

o"j being the Pauli matrices in spin space. The quantities Ao and Aj in 
Eq. (||) determine the singlet and triplet components of the order parameter 
A, respectively. In what follows, we shall assume that the structure of Ao 
and Aj is such that Ao(p) = Ao(p), Aj(p) = it^p^Afj?), where Rik is a 
real rotation matrix. The wave function of a Cooper pair for singlet-triplet 
pairing reads as 

g(p) = 9o(p)v2 + g{v)o<J2 (4) 

To derive the self-consistent equations it is necessary to specify the energy 
functional which we set in the form 

E(f,g) = E (f) + E int (g), E (f) = 2 £ e (p)/o(p), e (p) = (5) 

p 

Eint(g) = yY,9o(p) L s(p,p')go(p') + y^2g*(p)L t (p,p')g(p') 
pp' pp' 
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Here /o is the coefficient of the unit matrix do in expansion of the distribution 
function / in the Pauli matrices; L s ,Lt are the singlet and triplet anomalous 
FL interaction amplitudes with the structure L s (p, p') = L s (p,p'), Lt(p, p') = 
L t (p,p')p°p'° . The use of Eqs. ([[]), @,(||) allows to represent the self- 
consistent equations in the form 

a / \ l rr , JA±M + i. Y E + (q) , A_(g) ^ u Y Q E_(q) \ 

(6) 

K> 3 \^ + (g) 2 E_{q) 2 



Here A± = Ao ± A, = y£ 2 + A±. Eqs. (|6|) represent themselves the 
system of the integral equations for determining the singlet and triplet order 
parameters in the case of ST pairing of quasiparticles in a superfluid FL. 

Let us give an analysis of Eqs. (||), using the model representations of 
the BCS theory (the amplitudes L s and L t are not equal to zero only in a 
narrow layer near the Fermi surface: |£| < 6, 9 <C £f)- Then the quantities 
Ao = Ao(p = Pf), A = A(p = pf) can be found from the relations 

A = \x{l + d(x, T)), A = - d(x, T)) (7) 
where x is the solution of the equation 

d{x, T) ■ d(x • d(x, T),T) = D(x, T) = 1, (8) 

iM = WHcj , A(x , T)= /| tmh |, 

/ _ 0t v F L S)t {p = PF ,q= p F ) 



E = Je+x\ g> = ^, g Stt 



2 j ys,r — | 

One-gap solutions are obtained as solutions of the equations d(x, T) = 1 
(singlet pairing) and d(x,T) = — 1 (triplet pairing), while two-gap solutions 
correspond to the case d(x,T) ^ ±1 (singlet-triplet pairing). The critical 
temperatures of the singlet and triplet superfluid transitions are found from 
the equations d(0,T s ) = 1, d(0,Tt) = — 1 respectively. Assume for definite- 
ness, that g s > g' t . An analysis of the behaviour of the function D{x,T) 
shows (see Fig. ||) that equation D(x,T) = 1 has no one-gap solutions at 
temperatures T > T s , only one singlet solution exists for T < T < T s , while 
for T s t < T < T t the system is characterized by two one-gap (singlet and 
triplet) solutions. Finally, at T < T s t we have two new ST solutions in addi- 
tion to the previous solutions. For determining the critical temperature T s t 
of transition to the state with ST fermion pairing we have the equations 

d(x,T) = -l, xd' x (x,T) = 2 (9) 
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Fig. 1. Behaviour of the function D(x,T) at different temperatures. The 
critical temperatures of various phase transitions in the model case with g s = 
0.25, $f t = 0.2 and 6 = 0.01e F are as follows: T s = 0.0045e F , T t = 0.0033ep 
and T st = 0.0028e F . 



The first of these equations indicates that ST solutions are continuously 
branched off the one-gap triplet solution, while the second equation is the 
condition that the derivative D' x (x, T) vanishes at the branching point. Cal- 
culation of the second derivative D" x (x, T) in the critical point gives D" x (x s t, T st ) = 
0, i.e., the mechanism of branching of ST solutions is the formation of in- 
flection on the curve z = D(x,T st ) for x = x s t- 

The numerical analysis of the equations for determining Ao, A in the 
model case considered above is given in Fig. ||. Note that if the coupling 
constants satisfy the inequality g' t > g s then ST solutions branch off the 
singlet one-gap solution. 



3. SYSTEMS OF FERMIONS OF TWO SORTS 

The problem of studying of phase transitions in nuclear matter has en- 
couraged a great interest (see Refs. Bo and references therein). However, 
earlier the normal-to-superfluid phase transition in nuclear matter has been 
studied only. In principle, at further lowering of temperature a superfluid 
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Fig. 2. Order parameters Ao and A vs temperature. st(s+), st(s-) and 
st(t±) are notations for the dependences of the order parameters Ao and A 
in two pair (±Aq, A) of ST solutions of self-consistent equations. 



FL of nucleons can again become unstable passing into a new superfluid 
state. Thus, we are speaking of the superfluid-to-superfluid phase transi- 
tions in nuclear matter. The basic formalism is laid out in more detail in 
Ref. i, where the phase transitions to one-gap superfluid states of symmet- 
rical nuclear matter have been studied. The main aspect now is analysis 
of the phase transitions to multi-gap superfluid states. For simplicity we 
assume that the energy functional is invariant under rotations in the config- 
urational, spin and isospin spaces. Hence, superfluid phases (Cooper pairs) 
are classified by specifying the total spin of the pair, S = 0, 1, the isospin 
T = 0,1, their projections S z and T z on the z axis, and the orbital angular 
momentum L = 0, 1, 2, ... Each superfluid phase is described by its own set 
of order parameters: at S = 0, T = by the scalar order parameter Aoo 
(singlet-singlet pairing of nucleons), at5' = l,T = 0or5 = 0,T = lby the 
vector order parameter A^q or A fc, respectively (k = 1,2,3) (triplet-singlet 
and singlet-triplet pairing of nucleons) and at S = 1, T = 1 by the tensor 
order parameter Ajfc (i, k = 1,2,3) (triplet-triplet pairing of nucleons). 

Further we shall study the unitary states of superfluid nuclear matter, 
for which the product AA + is proportional to the product of the unit matri- 
ces in the spin and isospin spaces, AA + oc (JqTq. We consider as a particular 
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case two-gap unitary states, for which the order parameter reads as 

A(p) = A 3 o(p)o- 3 o-2T2 + A 03 (p)a 2 T 3 T2 (10) 

In this case the wave function of a Cooper pair describes the superposition of 
states with the triplet-singlet and singlet-triplet pairing of nucleons (TS-ST 
states) with the projections of total spin and isospin S z = T z = 0. 

To derive the self-consistent equations of a nucleon superfluid FL it is 
necessary to specify the energy functional, which we set in the form 

E(f,g) = E (f)+E int (g), Eo(f) = 4^T e (p)/oo(p), £o(p) = f^, (11) 

p 

2 

Eint(g) = 7^X!^3o(p)^i(p.q)53o(q) +5o3(p)^2(p,q)5o3(q)} 

p>q 

Here m is the effective nucleon mass, /oo is the coefficient of the product 
To<7o i n expansion of the distribution function / in the Pauli matrices, V±, V2 
are the anomalous FL interaction amplitudes, which have the symmetry 
properties Vi(— p, q) = V{ (p, — q) = Vi(p, q), i = 1,2. The quantity m 



contains account of the normal FL effects and represents itse 



f the mass of 



a free nucleon, renormalized by the normal FL interaction EL The use of 
Eqs. (IH) , (pi) , (lll|) allows to obtain the self-consistent equations in the form 



of Eqs. @, where A± = A30 ± A03, E± = y £ 2 + A^_ and it is necessary to 
do the substitution L s — ► V%, L t /3 — > V2. In what follows, we shall choose 
the Skyrme forces aa as the amplitudes of NN interaction 

Vi, 2 (p, q) = t (l ± x ) + U 3 n a (l ± z 3 ) + iti(l ± ^i)(p 2 + q 2 ), (12) 

where n is the density of symmetrical nuclear matter, ti, Xi, a are some phe- 
nomenological parameters. Note that the amplitudes V\,V2 contain no de- 
pendence from the parameters t2,X2, because the amplitudes V±, V2 are the 
even functions of the arguments p, q (see details in Ref. 11) . There are sets 
of parameters ti, Xi, a, which differ for various versions of the Skyrme forces 
(we shall use the Ska, SkM, SkM* and RATP potentials i as well as the SkP 
potential Q). 

The temperature dependence of the order parameters A30 = A 3 o(p = 
Pf), A03 = Ao3(p = pf) is determined from Eq. (0),® in which it is nec- 
essary to do the above mentioned substitution of the amplitudes as well as 
the substitution of the coupling constants g s — > g\ , g' t — > 52 • The results of 
numerical integration for the potential SkM* at n = 0.4fm -3 are given in 
Fig. |3|. It is seen that two-gap solutions branch off from one-gap singlet- 
triplet solution. The results of calculations for the potentials SkM, Ska have 
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Fig. 3. Order parameters A30, A03 vs temperature. tsst(ts+),tsst(ts-) and 
tsst(st±) are notations for the dependences of the order parameters in two 
pair (1LA30, A03) of TS-ST solutions of self-consistency equations. 



the form being analogous to the form in Fig. ||; for the potentials RATP and 
SkP the self-consistent equations have no two- gap solutions. 

Note that since the coupling constants depend on density of nuclear mat- 
ter, there exists the possibility of phase transitions in density to multi-gap 
states of superfiuid nuclear matter as well. It is demonstrated on Fig. ^| where 
the results of numerical determination of the order parameters A30 , A03 as 
functions of density at T = are presented. The results of calculations 
for the potentials SkM, Ska have the same form whereas for the potentials 
RATP, SkP the self-consistent equations have no two-gap solutions. 



4. CONCLUSION 

Thus, we have pointed out the possibility of phase transitions in super- 
fluid FL, consisting of fermions of one and two sorts, to two-gap superfiuid 
states, corresponding to the superposition of fermion pairings with different 
values of spin (and isotopic spin). The self-consistent equations for these 
states differ essentially from the equations of the BCS theory and contain 
the one-gap solutions as some particular cases. Phase transitions to multi- 



A.I. Akhiezer et al. 




gap superfluid states can take place either in temperature or density if the 
coupling constants in various pairing channels depend on n. Note, that 
mixed pairing of quasiparticles (but not in spin space) was considered in 
RefJ3 (superposition of d x 2_ y 2 and d xy states in orbital space for HTSC) 
and for superconductors with overlapping bands, first, in Ref.i and in some 
recent works (see, for example, Ref.0 and references therein). 



ACKNOWLEDGMENTS 

This research is supported by BMBF. 



REFERENCES 

1. A.I. Akhiezer, V.V. Krasil'nikov et al., Phys. Rep. 245, 1 (1994). 

2. A.I. Akhiezer, A. A. Isayev, S.V. Peletminsky et al., Zh. Eksp. Teor. Fiz. 112, 3 
(1997) [Sov. Phys. JETP 85, 1 (1997) ]. 

3. Th. Aim, G. Ropke, A. Sedrakian et al., Nucl. Phys. A 604, 491 (1996). 

4. M. Baldo, 0. Elgar0y, L. Engvik et al., Phys. Rev. C 58, 1921 (1998). 

5. M. Brack, C. Guet and H.-B. Hakansson, Phys. Rep. 123, 275 (1985). 



Theory of Multi-gap Superfluidity 



6. J. Dobaczcwski, I. Hamamoto et al., Phys. Rev. Lett. 72, 981 (1994). 

7. R.B. Laughlin, Phys. Rev. Lett. 80, 5188 (1998). 

8. H. Suhl, B.T. Matthias and L.R. Walker, Phys. Rev. Lett. 3, 552 (1959). 



9. D.F. Agterberg, V. Barzykin and L. Gor'kov, |cond-mat/9907331 



